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We present evidence for diffractive exclusive dijet production with an invariant dijet mass greater 
than 100 GeV in data collected with the DO experiment at the Fermilab Tevatron Collider. A 
discriminant based on calorimeter information is used to measure a significant number of events with 
little energy (typically less than 10 GeV) outside the dijet system, consistent with the diffractive 
exclusive dijet production topology. The probability for these events to be explained by other dijet 
production processes is 2 x 10 , corresponding to a 4.1 standard deviation significance. 

PACS numbers: 12.38.Qk, 24.10.Ht, 13.87.Ce, 12.40.Nn 



Hard diffraction was first observed about twenty years 
ago in the UA8 experiment at the CERN pp collider 
SPS [1] and has been studied extensively in several exper- 
iments: HI and ZEUS at the DESY ep Collider HERA @ 
and DO and CDF at the Fermilab Tevatron Q . At hadron 
colliders, hard diffractive events are identified by the sig- 
nature of a hard scatter in the presence of a region devoid 
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of any activity in the forward region of the detector or 
by tagging beam hadrons in the final state. Hard diffrac- 
tive events can be described by the exchange of a color- 
less object (Pomeron) 0]. Diffractively produced objects 
such as dijets, diphotons and \c charmonium can be ob- 
served in the detector together with Pomeron remnants. 
A subset of hard diffractive events in which both incom- 
ing hadrons remain intact is defined in such a way that 
all the energy not carried away by the outgoing beam 
particles is used to produce the diffractive system [3, 0] . 
This mechanism is defined as hard exclusive diffractive 
production. We search for this production mechanism in 
the sample of dijet events with large dijet invariant mass, 
corresponding to large values of the reduced center-of- 
mass energy of the Pomeron's system. 
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FIG. 1: Production of central dijet events in hard diffraction: 
(a) single diffraction, in which only either the proton or the 
antiproton is diffracted by a Pomeron (IP) exchange, while 
the other breaks up; (b) inclusive double Pomeron produc- 
tion, where proton and antiproton remain intact, and addi- 
tional QCD radiation can be observed from Pomeron rem- 
nants; and (c) exclusive diffractive production where both 
protons remain intact and only the dijet system is produced 
in the central region. 



Exclusive diffractive production (EDP) of a final state 
X, pp — >• p + X +p, has been proposed as a future search 
channel for new physics, as well as for the Higgs boson, 
at the Large Hadron Collider (LHC) Q- In this pro- 
cess, the kinematic properties such as the mass of the 
object X can be computed with high precision by mea- 
suring only the momentum loss of the outgoing protons 
in the final state. The CDF Collaboration reported the 
observation of exclusive diffractive events in the dijet, di- 
electron, diphoton and charmonium channels [8|. These 
results support the existence of EDP, but are restricted 
to low mass objects (typically less than 100 GeV), while 
at the LHC, searches for new physics are expected to ex- 
tend to higher masses. In this article, we report evidence 
for exclusive diffractive dijet production with invariant 
masses greater than 100 GeV in data collected by the DO 
experiment. 

We consider three different classes of hard diffractive 
production in addition to non-diffractive production: sin- 
gle diffractive (SD) dijet production (Fig.Q](a)), inclusive 
diffractive production through double Pomeron exchange 
(IDP) (Fig. [T](b)), an d exclusive diffractive dijet produc- 
tion (Fig. [I](c)). In SD, one of the beam hadrons remains 



intact while the other breaks up. In IDP, both beam 
hadrons are intact after the collision. The IDP and EDP 
processes with proton dissociations are expected to be 
suppressed by about a factor ten relative to the channel 
where the beam hadrons remain intact Q. The parton 
distributions of the Pomeron are taken from recent HI 
measurements [l(| and used to compute the diffractive 
jet production cross section at the Tevatron. An ad- 
ditional multiplicative factor (gap survival probability) 
11] of 0.1 is introduced to account for soft production 
of particles from the underlying pp events that populate 
the rapidity gaps 

The background to EDP in the dijet mass region con- 
sidered here originates from SD, IDP and non diffrac- 
tive (NDF) events which have either low multiplicity 
or small energy deposits in the forward calorimeters. 
Due to the steeply falling nature of these distributions, 
these backgrounds are expected to be small. NDF back- 
ground events are simulated using the pythia v6.202 12 1 
Monte Carlo (MC) generator with default settings and 
the diffractive (SD and IDP) backgrounds are determined 
using the POMWIG v2.0 [13] and FPMC vl.O @ generators 
respectively. An EDP of dijet events at the lowest order 
of QCD @ is simulated using fpmc through the exchange 
of two gluons. 

The data used in this analysis were collected with the 
DO detector in the period between August 2002 and April 
2006 at the Tevatron Collider at a center-of-mass energy 
y/s = 1.96 TeV. The DO detector is described in detail 
elsewhere [lij ]. For this analysis, the most relevant com- 
ponents are the central and forward calorimeters used 
for jet reconstruction and the identification of a rapidity 
gap devoid of any energy (above noise) in the calorime- 
ter, respectively. The DO liquid argon and uranium 
calorimeter is divided in three parts housed in indepen- 
dent cryostats covering the following regions in pseudo- 
rapidity: \r/\ <1.1 (central calorimeter), and 1.6 < |?7| < 
4.2 (two forward calorimeters) where ?y = — ln[tan(#/2)] 
and 6 is the polar angle with respect to the beam axis. 
Jets in EDP events are expected to be more central than 
in the other jet production processes, therefore both jets 
are required to be central with a rapidity \y\ < 0.8, where 
the rapidity is defined as y = 0.5 \n(E + p z )/(E — p z ) 
where E is the jet energy and p z is the momentum com- 
ponent of the jet along the beam axis. The forward region 
of the calorimeter is used to check for the presence of a 
rapidity gap on each side of the dijet system. 

The instantaneous luminosity used in this analysis is 
required to be in the range [5 — 100] x 10 30 enr 2 s _1 , 
where the contribution from two or more pp interactions 
in a single event is in general much less than 20%. This 
reduces the contamination of multiple interactions in the 
same bunch crossing to the rapidity gap selection. Data 
were collected using a inclusive jet trigger requiring at 
least one jet in an event to be above a pr threshold of 
45 GeV on the uncorrected energy, in order to to select 
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FIG. 2: Jet trigger efficiency as a function of the leading jet 
Pt (p J T )• For events with p^ > 100 GeV, the efficiency is 
close to 100% and no correction is needed. 



exclusive diffractive events in the region of dijet invari- 
ant mass above 100 GeV. Due to prescales imposed to 
avoid saturating the data acquisition system rate capa- 
bilities, the equivalent integrated luminosity of the sam- 
ple is about 30 pb _1 . By comparing the highest-pr jet 
spectrum with data collected with a trigger with a lower 
px threshold of 15 GeV, the trigger was found not to be 
fully efficient for jet px between 60 GeV and 100 GeV and 
the Monte Carlo events were reweighted with the trigger 
efficiency in this jet px range. The trigger efficiency as a 
function of jet px is shown in Fig. [2] 

Jets are reconstructed using an iterative midpoint cone 
algorithm [l5j with a cone size 1Z — y (At?) 2 + (A^) 2 = 
0.7, where tp is the azimuthal angle. The highest-pT and 
second-highest px jets are required to have px greater 
than 60 GeV and 40 GeV, respectively, and only dijet 
events with an invariant mass greater than 100 GeV are 
selected. To enhance the number of events without ad- 
ditional QCD radiation [6J, the two jets are required to 
be back-to-back in azimuthal angle tp, with a separation 
Atp > 3.1. A possible contribution of fake dijet events 
due to cosmic rays is suppressed by the requirement that 
the missing transverse momentum is less than 70% of the 
leading jet transverse momentum. 

The MC events are required to satisfy the same selec- 
tion criteria as the data. They are processed through 
a GEANT-based [l6| simulation of the DO detector re- 
sponse and the same reconstruction code as data. To 
simulate calorimeter noise and the effects of additional 
pp interactions, data events from random pp crossings 
are overlaid on the MC events, using data from the same 
time period as considered in the analysis. The MC events 
are weighted to obtain the same instantaneous luminosity 
profile as the data to have the same additional energy de- 
posits in the forward region of the calorimeter as in data. 
The sum of the number of NDF, SD and IDP events is 
normalized to data and their fractions are determined by 




FIG. 3: Dijet invariant mass distribution for MC and data. 
Good agreement between the MC simulation and data is 
found after applying jet energy scale corrections and scale 
factors corresponding to the trigger efficiencies, the luminos- 
ity profiles, and the MC normalization. 

the theoretical cross sections. The EDP contribution is 
expected to be negligible at this stage. Good agreement 
between the MC simulation and data is seen in the dijet 
invariant mass distributions in Fig. |3l By varying the re- 
quirement on the leading jet px, the uncertainty on the 
normalization was estimated to be 5%. 

To discriminate between exclusive events and back- 
ground (NDF, SD and IDP) we exploit the large rapid- 
ity gap that is expected between the central jets and 
the proton and antiproton beams. Two separate re- 
gions of pseudorapidity rj are defined in the calorime- 
ter far from the two central jets. The very forward re- 
gion (3.0 < \r/\ < 4.2) allows discrimination of diffractive 
events (SD and IDP) from NDF events, which are accom- 
panied by beam remnants in this region of the calorime- 
ter. The intermediate forward region (2.0 < \rj\ < 3.0) 
is used to identify EDP events, since they show larger 
rapidity gaps than SD and IDP. To prevent noise con- 
tamination in the calorimeter region under considera- 
tion, noisy cells in the forward region of the calorime- 
ters, which present an occupancy that differed by more 
than five standard deviations from the average, are re- 
moved. The cell response in MC was also adjusted to 
data by applying a MC-to-data correction factor for each 
cell. This correction factor was obtained using data col- 
lected requiring either minimal activity in the DO lumi- 
nosity counters or the presence of low px jets. After 
performing these corrections, the calorimeter cell infor- 
mation was used to form the following variable in order 
to discriminate between the different classes of events 

A=iexp(- ]T £ T ) + iexp(- £ E T ). (1) 

2<|r;|<3 3<|r)|<4.2 

Here Ex is the transverse energy in a given cell, and the 
sum is performed over all cells in the indicated rapidity 
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FIG. 4: Distribution of A normalized to unity for all MC 
samples. EDP peaks at A > 0.85. The EDP contribution at 
low A values is due to pile-up events, where a second proton- 
antiproton inelastic scattering occurs in the same bunch cross- 
ing. 



range. Figure [4] displays the A distribution normalized 
to unity for all MC samples. Also shown is the expected 
distribution from EDP events, showing a characteristic 
peak at A > 0.85, corresponding to energy deposits in 
the forward calorimeter which are typically smaller than 
10 GeV. 

Systematic uncertainties are assessed on the MC back- 
ground prediction of the differential distribution for 
A. The leading systematic uncertainty is due to the 
calorimeter cells calibration factors. They are varied 
simultaneously by three standard deviations from their 
central value leading to a change of 25% of the back- 
ground for A > 0.85. The effect of the jet energy scale 
uncertainties modifies the background by 12%. The jet 
energy resolution in simulation has been varied to match 
the data, yielding a small change of the normalization of 
0.5% which is assigned as an uncertainty. To estimate the 
uncertainties of the trigger efficiency correction and the 
instantaneous luminosity reweighting, the analysis was 
repeated using a 15 GeV jet pr trigger threshold result- 
ing in a 3% systematic shift. An additional systematic 
uncertainty due to the MC to data normalization is es- 
timated to be 5%. An uncertainty of 50% on the SD 
and IDP MC cross sections accounts for the uncertainty 
on the partonic structure of the Pomeron and survival 
probability gap factor. The non-diffractive parton distri- 
bution function uncertainties were considered and found 
to be negligible with respect to the other uncertainties. 
The total background prediction is 5.4 t.2'9 events and 
26 signal candidate events are observed in data. 

Figure [5] shows the comparison of the A distributions 
in data and MC (NDF, SD and IDP) normalized to 
their leading order cross sections. Good agreement is 
observed between data and MC except at high values of 
A where EDP dominates. The significance of the ex- 
cess with respect to the NDF, SD and IDP backgrounds 
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FIG. 5: (a) Distribution of A for data and contributions 
from NDF, SD and IDP. The hatched band indicates the total 
uncertainty on the background, (b) Normalized difference 
between data and NDF, SD and IDP processes. The red solid 
lines are ± 1 standard deviation systematic uncertainties on 
the background. 



is determined using a modified frequentist method [17]. 
It is obtained via fits of the signal+background and 
background-only hypotheses to pseudo-data samples con- 
taining only background. The effect of systematic uncer- 
tainties is constraintcd by maximizing a likelihood func- 
tion for background and signal+background hypotheses 
over all systematic uncertainties. Pseudo experiments 
used to determine the significance of the EDP signal in- 
clude variations over each systematic uncertainty. The 
observed significance corresponds to the fraction of out- 
comes that yield an EDP cross section at least as large 
as that measured in data. Seven bins are used as input 
for the significance calculation: six bins for A between 
0.1 and 0.85, where the predominant region used in the 
MC normalization is removed, and the A > 0.85 bin. 
The probability for the observed excess to be explained 
by an upward fluctuation of the background is 2 x 10~ 5 , 
corresponding to an excess of 4.1 standard deviations. 
Table U gives the observed number of events compared 
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FIG. 6: Dijet invariant mass distribution for stacked MC 
(NDF, SD and IDP) and data after applying the requirement 
on A > 0.85. The EDP distribution is shown without adding 
the stacked background. 



to background and EDP expectations. Figure [6] displays 
the dijet invariant mass distribution for A > 0.85. To 
illustrate the differences between the diffractive dijet ex- 
clusive events with A > 0.85, where the calorimeter has 
little energy deposition outside the central region, and 
the non-diffractive events, two event displays are shown 
in Fig. [7J The excess in data can contain events where 
the proton is dissociated into low-mass states that escape 
detection. The contribution of such events is estimated 
to be up to 10% of the EDP cross section 

To summarize, we have presented evidence at the 4.1 
standard deviation level for events consistent with the 
exclusive dijet production event topology in pp collisions 
at a center-of-mass energy y/s = 1.96 TeV at high di- 
jet invariant mass (Mjj > 100 GeV). These are the 
highest mass states studied for exclusive production in 
hadron colliders. Such event signatures are expected to 
play an important role in future studies at the Tevatron 
and LHC. 

We thank the staffs at Fermilab and collaborating 
institutions, and acknowledge support from the DOE 
and NSF (USA); CEA and CNRS/IN2P3 (France); 
FASI, Rosatom and RFBR (Russia); CNPq, FAPERJ, 
FAPESP and FUNDUNESP (Brazil); DAE and DST (In- 
dia); Colciencias (Colombia); CONACyT (Mexico); KRF 
and KOSEF (Korea); CONICET and UBACyT (Ar- 
gentina); FOM (The Netherlands); STFC and the Royal 
Society (United Kingdom); MSMT and GACR (Czech 
Republic); CRC Program and NSERC (Canada); BMBF 
and DFG (Germany); SFI (Ireland); The Swedish Re- 



[9 

[10 
[11 

[12; 

[13 
[14 

[is; 
[ie; 

[17 



UA8 Collaboration, R. Bonino et al, Phys. Lett. B 211 
(1988) 239. 

HI Collaboration, A. Aktas et al., Eur. Phys. J. C 48, 715 
(2006); ZEUS Collaboration, S. Chekanov et al., Nucl. 
Phys. B 831, 1 (2010); ZEUS Collaboration, S. Chekanov 
et al., Nucl. Phys. B 713 3 (2005). 

DO Collaboration, B. Abbott et al, Phys. Lett. B 531, 52 
(2002); DO Collaboration, B. Abbott et al., Phys. Lett. B 
574, 169 (2003); CDF Collaboration, D. E. Acosta et al., 
Phys. Rev. Lett. 91, 011802 (2003); CDF Collaboration, 
D. E. Acosta et al., Phys. Rev. Lett. 93, 141601 (2004). 
G. Ingelman, P. E. Schlein, Phys. Lett. B 152 (1985) 256; 
M. Boonekamp et al., Acta Phys. Polon. B 40 (2009) 2239 
and references therein. 

A. Schafer, O. Nachtmann and R. Schopf, Phys. Lett. 
B 249, 331 (1990); J. D. Bjorken, Phys. Rev. D 47, 101 
(1993). 

V. A. Khoze, A. D. Martin and M. G. Ryskin, Eur. Phys. 
J. C 14, 525 (2000); V. A. Khoze, A. D. Martin and 
M. G. Ryskin, Eur. Phys. J. C 19, 477 (2001) [Erratum- 
ibid. C 20, 599 (2001)]. 

M. Boonekamp, R. B. Peschanski and C. Royon, Phys. 
Lett. B 598 (2004) 243. 

CDF Collaboration, T. Aaltonen et al., Phys. Rev. D 77 
(2008) 052004; CDF Collaboration, T. Aaltonen et al., 
Phys. Rev. Lett. 99 (2007) 242002; CDF Collaboration, 
T. Aaltonen et al., Phys. Rev. Lett. 102 (2009) 222002; 
CDF Collaboration, T. Aaltonen et al., Phys. Rev. Lett. 
102 (2009) 242001; CDF Collaboration, A. Abulencia et 
al., Phys. Rev. Lett. 98 (2007) 112001. 
J. Monk and A. Pilkington, Comput. Phys. Commun. 
175, 232 (2006); M. Boonekamp et al., in Proceedings 
of the Workshop of the Implications of HERA for LHC 
physics, DESY-Proc-2009-02. 

HI Collaboration, A. Aktas et al., Eur. Phys. J. C 48 
(2006) 715. 

M. G. Ryskin, A. D. Martin and V. A. Khoze, 
Eur. Phys. J. C 60, 265 (2009); E. Gotsman et al., 
|arXiv:hep ph/0511060| 

T. Sjostrand et al., Comput. Phys. Commun. 135, 238 
(2001). 

B. E. Cox and J. R. Forshaw, Comput. Phys. Commun. 
144, 104 (2002). 

DO Collaboration, V.M. Abazov et al., Nucl. Instrum. 
Methods Phys. Res. A 565 (2006) 463. 
G.C. Blazey et al., in Proceedings of the Workshop: 
QCD and Weak Boson Physics in Run II, edited by 
U. Baur, R.K. Ellis, and D. Zeppenfeld, Fermilab-Pub- 
00/297 (2000). 

R. Brun and F. Carminati, CERN Program Library Long 

Writeup W5013, 1993 (unpublished). 

W. Fisher, FERMILAB-TM-2386-E (2006). 



8 



Sample 


NDF 


IDP 


SD 


EDP 


BKG 


DATA 


All A 


409527 


48.3 


2930 


30.9 


412505 


412505 


A > 0.85 


-2.9 


u - y -0.5 


o.2 ±»-.i 


12.9 ±J5 


5 4 +4 - 2 

-2.9 


26 



TABLE I: Number of predicted events for each MC sample for all A and for A > 0.85. The total uncertainties are quoted. 




(a) (b) 

FIG. 7: Event displays showing Et in the r\ — tp plane: (a) Exclusive diffractive event candidate: No energy deposition is 
present in the forward regions, only two central jets are observed in the detector, (b) Background event: In addition to the 
two jets present in the detector, energy deposition is present in the forward regions. The different colors correspond to energy 
deposits in different layers of the calorimeter. 



